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U ■ ABSTRACT 
< ■ 

We present mid-IR spectral decomposition of a sample of 48 Spitzer-selected 
ULIRGs spanning z~l-3 and likely Ljr ~ 10^^ -10^^ L©. Our study aims at 
quantifying the star-formation and AGN processes in these sources which recent 
1-^ ■ results suggest have evolved strongly between the observed epoch and today. To 

I ■ do this, we study the mid-IR contribution of PAH emission, continuum, and 

^ ■ extinction. About 3/4 of our sample are continuum- (i.e. AGN) dominated 

^ ■ sources, but ~60% of these show PAH emission, suggesting the presence of 

star-formation activity. These sources have redder mid-IR colors than typical 
optically-selected quasars. About 25% of our sample have strong PAH emission, 
but none are likely to be pure starbursts as reflected in their relatively high 5 /im 



> 

in 



^ ■ hot dust continua. However, their steep 30 /xm - to - 14 /xm slopes suggest that 

star-formation might dominate the total infrared luminosity. Six of our z ~ 2 
-H- . sources have i?PV6.2^0.3 /xm and Lu^m >W^hQ (implying Lir>10^3Lo). At 

these luminosities, such high EWq,2 ULIRGs do not exist in the local Universe. 
O ■ We find a median optical depth at 9.7/im of (rg.y^m) = 1-4. This is consistent 

with local JRAS-selected ULIRGs, but differs from early results on SCUBA- 
^ . selected 2; ~ 2 ULIRGs. Similar to local ULIRGs about 25% of our sample show 

extreme obscuration (rg.y^m >3) suggesting buried nuclei. In general, we find that 
our sources are similar to local ULIRGs, but are an order of magnitude more lu- 
minous. It is not clear whether our z ~ 2 ULIRGs are simply scaled-up versions 
of local ULIRGs, or subject to fundamentally different physical processes. 
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Introduction 



The pioneering work from the Infrared Space Observatory, ISO (see iGenzel fc Cesarsky 



2 pool , for a review) has already shown the power of mid-IR spectroscopy to reveal the domi- 
nant power sources among local infrared luminous galaxies. Diagnostic diagrams separating 



starburst from AGN doniinated sources have been developed from several studies (ILutz et al. 



19961. Il99a iGenzel et al.lll998l : iRigopoulou et al.lll999l : iLaurent et al.ll2000l : iTran et al 



2001 



Peeters et al.ll2004j ). Today, this work is being extended using the po werful Infrared Spec 



trograph (IRS; Houck et al. 2004) on the Spitzer Space Telescope (IWerner et al.l 120041 ). 



Its unprecedented sensitivity and efficiency have made it possible to expand the ISO work 
to much larger samples coveri 
higher spectral resolution (e.g 



Smith et al. 


2006; 


Brandl et al. 


2006; 


Schweitzer et al. 


2006 



Wu et al.ll2006l ). High S/N ratio IRS spectra from all these local samples have enabled var- 
ious teams to explore new param eter space and cr e ate new diagnostic tools characterizing 
different galaxy popula tions (e.g. 
20061 : ISpoon et al.lboOTh . 



Dale et al.l l2006l : lArmus et all l2006l . l2007al : ISturm et al. 



More importantly, IRS has enabled for the first time the mid-IR spectroscopy of IR 
luminous galaxies at cosmologically significant distances. Although limited to primarily low 
resolution data, the IRS spectra of high- 2; samples allow us for the first time to apply the spec- 
tral diagnosti c power of the mid-IR regi me to galaxies at 2; ~ 1 - 3, the epoch when LIRGs and 



ULIRGs (see ISanders fc Mirabellll996l . for a review) make a significant con tribution to the 



globa l averaged luminosity density and to the cosmic infrared background ( 



2005 



2005 



Le Floc'h et al. 



Dole et all l2006l : ICaputi et al.M2007n . Recent studies flHouck et al.l 120051 : lYan et al. 
20071 ) have produced such low resolution, mid-IR spectral samples of Spitzer-selected 
z~2 ULIRGs. Other IRS samples include near-I R-selected starbursts and X-ray selected 
AGN (IWeedman et al.ll2006l ) , and sub- mm galaxies (IMenendez-Delmestre et al.ll2006l : iPope et al 
2006bl ). Physical properties derived from these low- resolution, mid-IR spectra will allow us 
to put various disparate high-2 ULIRG samples in context, and test our understanding of 
IR luminous populations and their evolution. 

This paper is part two of a series of papers which present the results from a Spitzer 
spectroscopy survey of Spitzer-selected high-redshift galaxies (lYan et al.l l2005l ) . The first 
paper (PAPERI, Yan et al. 2007) describes the survey sample, data reduction, observed 
mid-IR spectra, and the basic properties of the galaxies in the sample, including redshift 
distribution, mid-IR luminosities and spectral types. This work (PAPERII) presents spectral 
decomposition, provides diagnostic tools, establishes the physical properties of our z ~ 2 
Spitzer selected ULIRGs, and discusses related limitations of current IRS spectra. The 
far-IR and radio properties of our sample are presented in Sajina et al. (2007, in prep). 
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Throughout this paper we adopt the cosmological model with f2M = 0.3, Q,a = 0-7, and 
Ho = 70. 



2. Data 

2.1. IRS spectral sample 

PAPERI gives the detailed description of the sample selection, presents the observed 
spectra and redshift measurements. The full sample has 52 sources selected to have 24 fim 
flux density brighter than 0.9mJy, and very red 24-to-8/im and 24-to-R color^. The IRS 
low resolution data revealed that 47 of the targets have measurable redshift from the features 
detected in their mid-IR spectra, with the majority (33/47 = 72%) at 1.5 < z <2 .7, and a 
smaller fraction (12/47 = 28%) at 0.65 <z< 1.5. An additional source (MIPS279) has a Keck 
redshift {z = 1.23) bringing our sample to 48. 



2.2. Broadband data 



We supplement the IRS spectra with Infr ared Array Camer a (IRAC; iFazio et al.l 120041 ). 



and Multiband Photometer for Spitzer (MIPS; Rieke et al. 20041 ) 70 /xm broadband photom- 
etry. 

The IRAC observations were obtained as part of the Spitzer Extragalactic First Look 
Survey (XFL^!!) dLacv et al.ll2005l ). The rms noise within a 6 " diameter aperture is 2.8, 3.2, 



15, 14.4yuJy at 3.6, 4.5, 5.8 and 8.0 //m respectively. We measured the aperture fluxes of our 
48 sources through apertures with a diameter of 7.3" (6 pixels) using the IRAC mosaiced 
images from Lacy et al. (2005). We applied aperture corrections of 1.112, 1.113, 1.125, and 
1.218 for the 3.6, 4.5, 5.8, and 8.0yum fluxes respectively. To obtain the sky variance, we 
used 12 sky apertures at > 10 pixels away from the source. From our sample, 5 sources are 
confused by nearby stars/bright galaxies. For these cases, we subtract the contaminant from 
the images by PSF fitting before performing the aperture photometry on our sources. The 
sources affected are: MIPS289, MIPS42, MIPS22661, MIPSllO and MIPS227. 



The MIPS 70 /im fluxes of our sample come from the combined scan map data (IFrayer et al 



20061 ) and pointed photometry data (PI: P. Choi). Most of our targets have roughly 700 sec- 



1(1) i?(24,8) = log,o{iyM2ifim)/,y,U{8pim) > 0.5; (2) i?(24,0.7) = logio(i//,(24/im)/z//,(0.7Mm) > 1.0. 
^ttp : //ssc . spitzer . caltech. edu/f Is/] 
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onds integration time per pixel. The 70 /im fluxes are estimated using a similar method 
as described for the IRAC bands. The aperture fluxes are measured through 3 pixel radius 
apertures. The total fluxes are computed by applying an aperture correction of 2.0. We also 
estimate flux errors using 12 sky apertures, and find typical rms flux errors of 1.3 mJy. All 
broad fluxes from 3.6-70/im are listed in Tabled! 



3. Fitting approach 



3.1. Spectral decomposition model 

The mid-IR spectra of a wide range of sources can be thought of as a superposi- 
tion of starburst-po\ yered PAH component and a predomiriantly AGN-power ed continuum 
component (see e.g. iLaurent et al.ll200d : ISturm et al.l l2000l : iTran et al.l 120011). This PAH- 



to-contiuum dia gnostic has been applied to ULIRG spectra as well (IGenzel et al.l Il998 



Tran et al.l l200ll ). However, lately the significant obscuration of mid-IR continua of the 
majority of ULIRGs has been highlighted as an ir nportant diagnostic parameter as it can 
dilute, or even erase, the power source signatures (ILevenson et al.l 120071 : ISpoon et al.ll2007l : 
Imanishill2007l ). With this in mind, below we describe our simple phenomenological approach 
aiming at distinguishing the contributions of PAH emission, continuum emission, and mid-IR 
obscuration. 

Our spectral fitting consists of two parts. The first is to determine continuum and 
opacity by fitting simultaneously three components: 1). an extinction-free PAH emission 
spectrum, 2). a power-law continuum, and 3). a mid-IR extinction law applied to the 
continuum. Because our spectra do not necessarily have the same relative PAH ratio as 
the assumed PAH template, the second step is to fit all available PAH features simultane- 
ously using Lorentz profiles, with the pre-determined continuum fit from the first step. The 
analytical formula describing the first fitting is as follows. 



Fu = apAHF„^PAR + aci^ "e (1) 

where opah and Oc are the amplitudes of the two components, F^ pp^n is a PAH template, 
a is the spectral slope, and t,^ is the extinction curve. The normalized extinction curve is 
fixed, with rg.r being the free parameter. We perform a simultaneous fit on the four free 
parameters above. 

We limit the fitting to the restframe 5-15/im spectral window which is the typical 
restframe coverage of our IRS spectra. The power-law approximation is assumed to be valid 
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in this limited region, but is expected to break down at shorter and longer wa velengths due 
to dus t sublimation and influence of colder dust components respectively (e.g. 
2007ah . 



Armus et al. 



Due to the range of redshifts (0.6 < z <3.2) of our sources, we see intrinsically different 
parts of the PAH spectrum. However, for comparisons within our sample, we need consistent 
continuum estimates regardless of the rest-frame spectral coverage. This requires adopting 
a single PAH template, and ignoring PAH variability at this point (see above and § 13. 2p . 
The pure PAH sp^trung (-PpAH)..is de rived from the IRS, low resolution spectrum of the 



starburst NGC7714 (IBrandl et al.ll2006l ). where we have subtracted its continuum as shown 
in Figure [1] left. 



We tested the effect of our choice of PAH template on the fit results (see § 14.11) . We 
looked at variations in the when substituting the NGC7714-based template with the 
M82 spectrum wh ich intrinsically ha s much less ~ 5 - 15 /im continuum emission than the 
average starburst (IBrandl et al.ll2006l ). and therefore we took it 'as is' without subtracting a 
continuum. We found the variations in the to be negligible, i.e. well within the continuum 
uncertainties. Therefore, for our data, the PAH-continuum decomposition is robust against 
different PAH template choices. 



•^This is referred to as a PAH spectrum simply because of the lack of continuum. It is more accurately 
thought of as "continuum-subtracted starburst spectrum". For simplicity however, we will refer to it as a 
"PAH template". 
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0.10 




Fig. 1. — Left: The solid curve is the IRS spectrum of NGC7714. The dashed hne is the a = 3 
power law that best fits the ~5-15/im continuum without extinction. The PAH template 
used is the residual (dotted hne). Right: The Milky Way ayerage (i.e. disk) extinction curve 
and the Galactic Center extinction curve (jChiar &: Tielensll2006l ). Here we normalized both 
to the K-band to highlight their differences in the ~5-15 fim regime of interest. 
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We choose the Galactic Center (GC) extinction curve (jChiar &: TielensI l2006l ) as hkely 
closer to the extreme environments in ULIRGs than a standard ISM extinction curve. Fig- 
ure ^ight s hows the GC extinctio n curve compared with the Milky Way disk (MW) extinc- 
tion curve Jchiar fc TielensI booeh . 

We assume screen extinctioE0, which means we ignore radiative transfer effects: absorp- 
tion and emission along the line of sight in an optically thick medium. 



This approa ch is similar to some previous studies of local ULIRGs (ITran et al.l 12001 



Lutz et al.ll2003[ ). 



3.2. PAH profile fitting 

To accurately measure the strength of observed PAH features, we simultaneously fit 
a set of Lorentz profiles of the continuum subtracted spectrum at approximately 6.2, 7.7, 
8.6, and 11.3/im|j. An example is presented in Figure [2j Below we give the expressions 
for line flux and equivalent width, as well as their associated uncertainties. The amplitudes, 
widths, and exact central wavelengths of the features are all free parameters. T he widths are 



not f ixed because of intrinsic variability such as due to secondary features (e.g. ISmith et al. 



20061 ). The central wavelengths are allowed to vary within a limited range largely due to the 
redshift uncertainties of our sample (see PAPERI). Denoting the flux density as F\ and the 
continuum flux density as Fc, the total flux and its associated uncertainty is determined by 
Equation [2l Note that these are the flux densities of each fitted component rather than of 
the original spectrum. 



* Since this is the most cfBcicnt way to absorb hght (compared with mixed and/or clumpy distributions), 
the optical depths derived in this manner are lower limits. 

^We ignore the 12.6 /im feature since it is blended with the Neii 12.8 /im line. 
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Fig. 2. — An example of our PAH fitting procedure. The best-fit continuum is shown in 
blue, the individual Lorentzian fits are shown as red dashed lines, while the solid red line is 
the total. For comparison, the dotted curve shows a cubic spline continuum, which misses 
most of the 7.7/im feature as determined using our continuum model (see also. ISmith et al. 
20061 ). The bottom panel shows the residual. 
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(2) 

where the hmits, Ai and A2 are such that they encompass 95% of the hne profil^. 

Assuming a constant rms and ignoring the continuum uncertainty, the above formula 
simphfies to the more famihar: 1 a = ^/iVpi^rmsAA. Equivalent widths (EW) and their 
uncertainties are determined by Equation [31 

EW = Ell ^dX, 

(3) 



-X2 



dX 



2 



In cases where the spectral coverage includes a specific line but it is too weak for profile 
fitting, we estimate 1 a upper limits by using the uncertainty equations for the flux and 
equivalent width respectively. This requires adopting a minimum width, 7mm- For 6.2, 7.7 
and 11.3 lines respectively this is: 0.1 fim, 0.2 /im and 0.1 /xm. These values are consistent 
with the minimum widths found for the detected sources. 

Both the flux density and equivalent width uncertainties are affected by not just the rms 
of the spectra, but the uncertainty of the continuum level as well (see Equations [21 As 
we discuss in § 4.1, we use Markov Chain Monte Carlo (MCMC) fitting, which provides us 
with the probability distributions of the fit parameters. These include both the uncertainty 
in the data and the uncertainty purely due to the continuum level. We can derive the 
latter by subtracting the variance of the data from the the MCMC-derived variance by using 
a^gf. = a'^ (Fc) + a'^ (Fx) . This is necessary in order to avoid counting the data rms twicqj. 
This is a non-negligible effect as we find that the continuum error dominates at 7. 7 /im and 
11.3 fim , and for ~ 1/2 of the sources at 6.2 /xm . 



^For a Lorentz profile this is 57, where 1 7 is the half-width encompassing half the power in the line. For 
comparison, in a Gauss profile, this corresponds to 2 a. The Lorentz width parameter, 7, translates into the 
more famihar Gaussian FWHM via 7 FWHM/1.8 

^Note that this therefore is different from the uncertainty in the continuum luminosities, which we estimate 
from the marginalized probability distribution as is, since in that case we want to include the data rms as 
well (see S OD . 




Fig. 3. — Left: The IRS spectra (± la) with the best-fit continuum (blue) and contin- 
uum+PAH (red). Right: Including the IRAC and MIPS data (squares). Overlaid is our 
model plus the host galaxy (dashed green) and cold dust excess. 
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Fie;. 3. — Continued. 
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Fig. 3. — Continued. 
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4. Fitting results 

4.1. Spectral Decomposition 

We fit the model described in § l3.1l using Markov Chain Monte Carlo (MCMC; see Sajina 
et al. 2006 for a description and references). The left-hand panels of Fig. [3] show the 5- 
15 /im spectra (over which the fit was performed) overlaid with the best-fit model. For the 
z >2.5 sources, we include the MIPS70/xm point in the fit to help constrain the A> 13 /im 
continuum. A number of typically weaker feature^ are likely to be present in our data, 
but are unaccounted for by the model. These include the neon lines, especially at 12.8 /xm; 
absorption feature s at 6 .0 /im, 6.9 /zm and 7.3 /im due to water ice a nd hydrocarbons (se e 



Spoon et al.l |2002| . 120031 ) . as well as the 9.7 /im H2 emission feature (IHigdon et al.ll2006l ). 
With this caveat, we find that our simple model provides reasonable fits to the spectra (the 
fit for M1PS464 was forced to go through the 5.5 /xm continuum due to the strong 6 /xm 
absorption - see appendix). The best-fit continuum parameters are given in Table 1. The 
uncertainties are determined from the probability distributions sampled by the MCMC. 

Our model fitting is performed on the ~ 5 - 15 /im spectra alone as stated above. How- 
ever, the availability of the IRAC and MIPS70 /im data allow us access to additional diag- 
nostics as described in § 14.31 Therefore in the right-hand panel of Figure [3], we show the 
IRS spectra along with the available IRAC3.6/im, 4.5 /xm, 5.8 /im, and 8.0 /xm as well as 
MIPS70 /xm broadband points. We overlay the best-fit model shown in the left-hand panel, 
but with one important modification. Since the power-law approximation cannot physi- 
cally extend to much shorter and longer wavelengths, we truncate it below 5 /xm and above 
15 /xm§. It is expected that the IRAC points would be sensitive to host galaxy stellar light 
(if present), while the 70 /xm points (Arest^20 /xm) would be sensitive to cooler dust emission 
(if present). The data do not allow sophisticated study of either of these components, and 
hence we characterize each with a single template (see § 14.31 for details). For each of these 
two components a simple linear fit is performed to find the best-fit amplitude for the corre- 
sponding template. In the right-hand panels of Figure [3l we show the best-fit host galaxy 
component (dashed green curve), and the best-fit thermal dust component (dotted green 
curve). The red curve shows the best-fit mid-IR model with these two components added to 
it (see § 14.31 for a discussion). 



^The start and end couple of microns of the spectra frequently show apparently strong, but likely spurious 
features. 

^This is done by assuming the power-law is a collection of blackbodies and setting Tmin — 600 K and 
Tmax = 1000- 1500 K. 
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4.2. Mid-IR continuum slope and opacity 

Figure H] shows the distribution of optical depth vs. continuum slope for our sample. 
There is a notable lack of sources with low r and flat spectra (i.e. a ~ 1). This is due to one 
of our selection cuts, log(S'24/5'8) > 0.5 (see PAPERI). The effect of this cut as a function of 
optical depth is shown in Figure H] for z = 2 and 2; = 2.5. The median parameter values are: 
(a) =2.0 and (tq.t) =1.4. 
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Fig. 4. — The continuum parameters a and r. The curves show the effect of our S24/ Si 
color selection for 2; = 2 and z — 2.5 as labeled. 
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Our optical depth measurements assume an extinction curve acting on a power-law 
continuum. Another common measure is the apparent depth of the si licate feature com- 
pared with the co ntinuum defined by the 5/im and 15 fim points (e.g. iBrandl et al.l 12006 



Spoon et al.l 120071 ). Given a fixed extinction curve, the two quantities have a simple scaling 



relation. In the case of the Galactic Center extinction curve we adopt, this relation is simply: 
T9.7 = 1-4 X r^i obs- 



4.3. Monochromatic luminosities 

We use the best-fit continua of our IRS spectra to derive monochromatic luminosities 
(z/Lj,) in two pure continuum regions. We pick 5.8 /xm due to its sampling the hot dust 
continuum just shortward of the 6.2 fim PAH feature, and 14 /xm for its sampling the con- 
tinuum between the 9.7 /xm and 18 /xm SiO absorption features (it is a minimum in the GC 
extinction curve see Figured]). We read these luminosities from the best-fit continuum fits 
not directly from the data. This is necessary because of variable rest-frame spectral coverage, 
noise spikes at the beginning and end of spectra. This also means that observed absorption 
features (e.g. MIPS464) not present in the model, do not affect these luminosities. 

The IRAC and MIPS70 fim photometric points allow us to sample the SEDs of our 
sources outside of the regions covered by the IRS spectra (see § 14. ip . Here we quantify the 
emission shortward and longward of the available spectra by deriving the monochromatic 
1.6 /im and 30 /xm emission wherever possible. As discussed in § 14. H the IRAC points (rest- 
frame ~ 1 - 3 /xm ) are likely to be sensitive to direct stellar emission from the host galaxy, 
which has a characteristic bump at 1.6 /xm. We only accept that we are likely seeing host 
galaxy light if: 1) the source is detected in both 3.6 /xm and 4.5 /xm, and 2) if the spec- 
tral shape changes color from the dust continuum (e.g. yes for MIPS15949, but no for 
MIPS15958). These data do not allow us to explore the properties of the stellar popula- 
tions. All we can address is the amplitude of the host galaxy emission whenever present. 
To do so we subtract from the IRAC points the mid-IR continuum (truncated below ~ 5 fim 
as describ ed in ^14. ip. and perform a simpl e linear fit using a 10 Gyr-old stellar population 



template (IFioc fc Rocca-Volmerangd 119971 ). The choice of template here is irrelevant, as 
it merely allows us to find the amplitude of the host galaxy emission (quantified via the 
monochromatic 1.6 /xm luminosity). 

On the other hand ai z < 1.6, the MIPS70 /im point samples restframe ~30/xm. If it 
is present, this traces the colder dust emission dominating the typical star-forming galaxy's 
SED. We therefore fit a 50K thermal emission to the MIPS70/im point, whenever it is in 
excess of the truncated power-law continuum. We emphasize that, as in the case of the stellar 
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population, the current data do not allow us to determine the SED peak (and hence dust 
temperature). Therefore, the 50 K blackbody is used only as a spectral templates allowing 
us to derive the monochromatic 30 fim luminosity. 

The broadband-derived 1.6 /im and 30 /xm luminosities, along with the 5.8 /im and 
14 /im luminosities derived from the IRS spectral fitting itself are shown in Table [2l 



4.4. PAH features strength 

To determine the PAH features strength, we start by subtracting the best-fit continuum 
from the spectra. We then perform the simultaneous Lorentz profile fits to the main features 
as described in § 13. 2[ 

The PAH equivalent widths and fluxes are tabulated in Table [3l where the values less 
than 2a (this includes both the data and continuum uncertainties) are shown as upper limits 
(sources without spectral coverage have no value). We leave off the 8.6 /im feature as it is 
too confused with the 7.7/im feature and affected by the continuum uncertainties (due to 
the Si feature). We only include it in the fitting primarily for the sake of improving the 
7.7/im feature strength estimates. 

Out of a total of 48 sources, 38 have data at 6.2 fim of which 2 have the 6.2 /im PAH 
feature detected at >3(t (9 at >2(t). Out of 46 sources with spectral coverage, 24 have 
their 7.7/im feature detected at >3(t (36 at >2cr). For the 11.3 /im feature, out of 37 
sources with data 12 are detected at >3a (20 at >2cr). Although the uncertainties vary 
within the sample, the minimum detectable equivalent widths are roughly EW6.2^0.2/im, 
EW7.7 >0.3 /im , and EW11.3 >0.3/im. Care must be taken when comparing our EW mea- 
surements with those obtained using different approaches. The difference in 6.2 /im feature 
equivalent width measurements are well within ~ 50% between different fitting approaches. 



using a cubic spline continuum (e.g. 



Brandl et al. 


2006; 


Armus et al. 


2007a) 



The luminosity of the strongest PAH feature, L7.7 is given in Table [21 along with the 
key monochromatic luminosities tracing the SEDs discussed in § 14.31 
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Fig. 5. — The stacked spectra of weak-PAH sources divided into optically-thin (r9.7<l) 
and optically-thick {rgj > 1). Note that a simple model with our PAH template added to a 
continuum (in blue) is a reasonable fit to the data. The power-law slope is 2.1 for the low-r 
spectrum, and 1.7 for the high-r spectrum. A r9.7 = 2.0 extinction is applied to the later. 
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4.4- i- Stacking weak-PAH spectra 

In Figure[5], we show the stacked spectra of the weakest-PAH sources (i.e. EW7.7 < 0.7/im ), 
separated into optically-thin (rg.r < 1) and optically-thick (rg.y > 1) sources. Individual spec- 
tra were first scaled to have the same total power (i.e. ^ i^F^,) which is less biased than 
normalizing at a particular wavelength. The stacked spectra are the weighted means of these 
scaled spectra. These stacked spectra clearly show that the PAH+continuum interpretation 
is valid even if this is not always obvious in individual spectra. As discussed in § 14.11 other 
features are likely to be present in the spectra, which could lead to misidentification of 
features and therefore wrong redshifts. Figure suggests that this is unlikely for the vast 
majority of the sample. 

The stacked spectra of the high-r, weak-PAH sources highlights that the distorted hump 
at ~ 8 fim is easily explained by a standard PAH emission spectrum superimposed on a 
deeply obscured continuum. The presence of the 6.2 fim feature confirms this. However, we 
have excluded MIPS429 and MIPS464 due to their exceptionally strong ~ 6 /xm absorption 
features likely due to water ice and/or hydrocarbons (see Appendix A). The average of the 
remaining obscured weak-PAH sources still shows a weaker 6.2 /im feature compared with 
our template, which suggest that the same absorption might be present in these sources as 
well. On the other hand, the weak-PAH low-r sources also show a 6.2 /xm PAH feature 
that is weaker relative to our PAH template. This may indicate smearing due to redshift 
uncertainties, PAH obscuration, or variations in the intrinsic PAH features strength. 



5. Starburst-AGN diagnostics 
5.1. The PAH features strength 

5.1.1. 7.7 PAH feature 

We focus on the 7.7/im feature as it is usually the strongest and this spectral region is 
observed for all of our sample. 

The downside is that it is the feature most affected by the continuum model. The 
6.2 /im feature is much cleaner in that respect, a nd therefore is usually used as proxy fo r the 



PAH strength of local LIRGs and ULIRGs (e.g. lArmus et al.ll2007al : lOesai et al.l 120061 ). In 
Figure [6] we show the restframe 6.2 /xm equivalent width [EWq,2) vs. the restframe 7.7/im 
equivalent width {EWj^f). We find that i?PV6.2^0.2/im is the detectab i lity thr eshold for 
the 6.2 /im feature, which translates into i^PVy.y >0.8/im. lArmus et al. J2007ah find that 
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EWq,2 ^ 0.4 fim separates pure starbursts from comp osite systems. Although this is uncertain 



to ~50% due to the different fitting approach in lArmus et ahl (l2007al ). it suggests that 
our borderhne detected sources are hkely to include significant AGN contribution as well. 
Therefore, this equivalent width cut should not automatically be translated into a starburst- 
AGN separation, but merely into a separation into strong-PAH (largely corresponding to 
detected 6.2 /im feature when the data are available) and weak- PAH sources. From now on 
we refer to sources with E^Wy.y > 0.8 /xm as PAH-rich, and those with EW^/i <0.8 as 
PAH-poor. We emphasize that PAH-poor does not imply no-PAH since 67% of the PAH- 
poor sources also have detectable 7.7/im PAH feature, and likely have some star-formation 
activity in addition to the dominant AGN. A quarter of our sample are PAH-rich sources, 
and three quarters are PAH-poor sources. In the rest of this section, we attempt to attach 
more physical meaning to this empirical classification. In particular, we ask whether, or 
not, our strong-PAH sources are composite systems with comparable AGN and starburst 
contribution to the mid-IR light, or some of them are clearly starburst-dominated. 

Lastly, for the purpose of comparison, we have applied our fi tting approach on a small 
sub-sample of the local ULIRG sample from lArmus et aP (2007a). These are overplotted in 
Figure El and will be included whenever possible in subsequent figures. 



- 30 - 




tRAS^57^5®5 

^Jf^,K;i,ni,zL , , , 



I I I I I 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

EW6.2„,.[nm] 



Fig. 6. — Here we show the sources with EW7.7/im > 3 cr, where data are available 
at 6.2 /xm as well. The 6.2 /im detection limit is EW6.2/im >0.2/im , corresponding to 
EW7.7/im >0.8/im. Both limits are shown as dashed lines. We will use the 7.7/im equiv- 
alent widths in the rest of the paper as it is available for most of our sources. From now 
on, red crosses denote PAH-rich sources (EW7.7/im >0.8/im), while blue circles denote 
PAH-poor sources (EW7.7/im <0.8/im). 
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5.1.2. Trends with redshift and luminosity 



In PAPERI, we presented the redshift distribution of our sample, showing a z^2.1 
peak and a smaller z~0.8 peak. Given our 24 /xm flux selection, we expect the lower-2 
galaxies to have 10^^ Lq (typical local ULIRGs), while the higher-^ galaxies to have up to 
LjR ~ 10^^ Lq (potential HyLIRGs). Locally the AGN-fra ction in local UL IRGs rises with 
total infrared luminosity, dominating above ~ 3 x 10^^ L0 (ITran et al.ll200ll : ?). If this holds 
for our sample, we would expect the strong PAH sources to be among the lower-z sub-sample, 
while the higher- 2;, higher-LjR sources to be continuum-dominated. 



- 32 - 




(^EW7.7<0.8^m 



0.1 1.0 
EW7.7 [|im] 

Fig. 7. — For the z>1.5, we find a trend of decreasing EWt,-j with increasing L14. As 
usual, we mark the strong- PAH sources (i.e. EW7.7> 0.8/im) with red crosses, while ISiS^'^ 
below this limit are shown with open blue circles. All z < 1.5 sources (larger symbols) are 
significantly below the trend found for the higher- 2; sources. For the whole sample, 25% are 
strong PAH sources - for the low-2 {z < 1.5) sources this fraction is 42%, while for the high-^ 
sources it is only 23%. 
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Figure [71 shows that within a given redshift bin, there is a trend of increasing L14 with 
decreasing EWt,-^. This might be a trivial result as increasing mid-IR continuum by definition 
decreases the equivalent width. We need total infra red luminosities t o study trends in PAH 
strength with luminosity as seen for local ULIRGs (ITran et al.ll200ll ). 



1.0 




NGC6240 (Armus et al. 2006a) 

Average starburst from Brandl et al. (2006) 

^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ 

4 6 8 10 12 14 

?l [jim] 

Fig. 8. — The average spectra of four of our 2; ~ 2 strong- PAH sources, normalized at 7.7/im . 
Our sources appear very similar to the local ULIRG NGC6240. Both are similar to the 
starburst template, except for the excess power at A < 6 /xm . 
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A more intriguing aspect of Figure [7| is that although it is true that strong PAH 
sources are predominantly low-z, there is a population of equally strong PAH sources at 
Figure [S] shows the average of strongest-PAH {EWj,t>2 fim) of these z~2 sources 
(MIPS289, MIPS8493, MIPS16144, and MIPS22530), normalized at 7.7/im. These sources 
have Li4 ~ 10^ ^ L(d, implying Ltr ^ 10^^ Lq. For comparison we overlay the average starburst 
template from lBrandl et alj J200e ). as well as the local ULIRG NGC6240 (Lm ~ 7 x 10^° Lq, 
from our own fit) (lArmus et al.l l2006l ). Our z^2 strong- PAH sources have spectra very 



similar to NGC6240; however, our sources are about an order of magnitude more luminous. 
Lastly, both our average z ~ 2 spectrum and NGC6240 show excess emission at A < 6 fim 
relative to the average starburst. We return to this point in § 15.21 



5.2. Mid-IR colors: hot and cold dust 

The inner dust tori around AGN can be heated to near sublimation temperatures 
(~ 1000- 1500 K), which dominates the mid-IR spectra provided these hot regions are not 
obscured (e.g. Pier & Krolik 1992). Such emission peaks in the ~3-6/im region leading 
to (Type 1) AGN to have flatter ~5-15/im spectral slopes (a~l, e.g. 3c2 73) than the 



typica l starburst (a~3, e.g. NGC7714). The typical slopes of local ULIRGs (lArmus et al. 



2007al ) lie roughly between these values, and are consistent with our median slope of a~2. 
We find no difference in the ~ 5 - 15 fim slopes distribution of the PAH-rich and PAH-poor 
sources. 

The longer wavelength SED is dominated by cold dust of < 50 K associated with star- 
formation. In § 14.31 and Figure [31 we show that many of our sources show very warm spectra 
peaking at ~ 30 /xm , while others show an excess over the mid-IR continuum, consistent 
with the presence of cooler dust. These were used to derive z/Li,(30/im) (see Table [2]) 
for sources with z<1.6 such that the MIPS70/im point is not too redshifted away from 
restframe 30 /xm. Figure [9], shows that while the ~ 5 /im-to-15/im colors are the same for 
both PAH-rich and PAH-poor sources (as is the case for the entire sample), the ~15/im- 
to-30 fim colors separate the two groups very well. Our PAH-poor sources are redder than 
the optically-selected quasars consistent with their being dustier. Our PAH-rich sources are 
generally consistent with the starbursts' colors although they show slightly bluer ~ 5 - 15 fim 
slopes suggesting enhanced hot dust. 
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Fig. 9. — The mid-IR color-color diagram, where only z<1.6 sources are shown as reliable 
30 fim luminosities cannot be obtained for the rest. For comparison we overlay local samples 
of starbursts and optically-selected quasars. The dashed line is to guide the eye. Note that 
among the reddest in both colors sources, the two populations are indistinguishable. 
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Fig. 10. — The z/L,^(5.8) vs. Lpauj.y- For comparison we overlay local samples of optically- 
selected AGN, normal star-forming galaxies, and starbursts. The exact positions of these 
are uncertain as weM4"^e converted the published PAH 7.7 /im feature strength values to our 
continuum definition. We have fitted the named sources with our method. 
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We expand on the last point, now using the entire sample, in Figure [TOl where we show 
the hot dust continuum luminosity z/L,^(5.8) vs. the 7.7/im PAH luminosity. For lower 
luminosity star-forming sources, the two are well correlated with slope unity suggesting at 
least at moderate luminosities, the mid-IR SEDs of purely star-forming galaxies are fairly 
constant. As expected, optically-selected quasars have stronger 5.8 /im continuum per 7.7/im 
luminosity. Our weak-PAH sources, at higher luminosities, are consistent with the local 
AGNs trend, despite their redder slopes. Although, to a lesser extend, our strong-PAH 
sources are also noticeably shifted with respect to the local starbursts. There are two possible 
explanations for this: 

1) They have significant AGN contributions leading to enhanced 5.8 fim continuum. 

2) The local starburst trend does not extend to the highest luminosities regime, but 
curves toward our strong-PAH sources making them starburst-dominated as well. In other 
words, extreme starbursts generate enhanced hot dust emission compared with lower lumi- 
nosity ones. 

The first explanation is very likely given the extreme luminosities of our sources. As dis- 
cussed in § 15.1.11 our equivalent width cut is more an upper limit on the starburst-dominated 
sources, since especially borderline sources are likely to include significant AGN contribu- 
tion as well. Therefore, we also look individually at the z ~ 2 sources with strongest PAH 
emission {EWj,-j > 2 /im , or EWq,2 ^ 0.3 /xm ). Even these strongest PAH sources, show excess 
emission at A<6/im, relative to the typical starburst (see Figure [H]), indicating higher hot 
dust contribution. Figure [8] also shows that our s trong-PAH sources are similar to NGC6240 
in that respect, which was already interpreted by lArmus et al.l (120061 ) as having ~ 20% AG N 
contribution to the bolometric luminosity (~ 50% of the mid-IR emission; iLutz et al.l (120031 )). 



On the other hand, hot dust emission has been observed to be associated with star- 



formation activity as well (iHunt et al.ll2002l : iLu et al.ll2003l ). An extreme exan aple, are high 
SFR intensity dwarf galaxies which show continuum-dominated mid-IR spectra (iGalliano et al 
20031 ) . Therefore, the possibility of enhanced hot dust emission being associated with ex- 
treme starbursts cannot be excluded. 



We conservatively adopt the first explanation above, but stress that the properties of 
high-luminosity, pure-starbursts remain unclear. 
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5.3. Near-IR host galaxy light 



It has been proposed that the presence of the 1.6 /xm bump, characteristic of stellar 
populations older than ^ 100 Myrs, is an efficient way of separat ing starburst- from AGN- 



dominated sources (e.g. IWeedman et al.ll2006l : iTeplitz et al.l 120071 ). The rationale is that the 



hot dust continuum of AGN overwhelmes their host galaxy's stellar bump. 
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vL„(1.6nm) [LqI V L„(1.6nm)/v L„(5.8nm) 

Fig. 11. — Left: The 1.6 /xm stellar bump strength vs. the 7.7/im PAH strength. The 
normal starforming galaxies are used to define the best-fit line shown. Note that both our 
PAH-rich and PAH-poor sources agree very well with this relation. Optically-selected AGN 
light to the right of this, suggesting that for the bulk of our PAH-poor sources, the 1.6/im 
emission is more likely to be due to host galaxy emission rather than that of the AGN itself 
as in the PG QSOs sample. Right: The 1.6/im stellar bump strength vs. the 7.7/im PAH 
strength normalized to the 5.8 /im luminosities. 
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Figure [3] shows that the IRAC points (which probe the 1.6 /im stellar bump at 2;~2) 
often (in 73% of the sources) have excess emission over the extrapolated mid-IR continuum. 
This excess is most likely the 1.6 /xm stellar bump of old stellar populations together with 
some unknown young stars contribution (see § 14.31) . The left-hand panel of Figure [TT] shows 
the 1.6 /xm luminosity vs. the PAH luminosity. Both our PAH-rich and PAH-poor sources 
(where the excess is observed) obey the relation defined by the local star-forming galax- 
ies. This suggests that even in the continuum-dominated sources (with a few exceptions) 
we do not see the AGN itself (accretion disk emission), but the host galaxy light alonq^. 
Circumstantial support is provided by optically-selected quasars whose 1.6 /im luminosities 
(McLeod & Rieke 1994a,b) lie significantly to the right of the star-forming galaxies relation. 
Although not shown for clarity, the AGN-subtracted H-band luminosities of the PG QSOs 
(McLeod & Rieke 1994a,b) agree much better with the star-forming galaxies relation sug- 
gesting the offset is predominantly due to the effect of AGN-dominance at 1.6 fim rather 
than suppressed 7.7 /im emission. 

The last point is supported by the right-hand panel of Figure [TT] shows z/Li,(1.6 /im ) 
vs. the 7.7 fim PAH feature luminosity both normalized to the 5.8 yum continuum. The 
effect of increasing hot dust emission is now obvious, as there is a clear trend between the 
starforming galaxies, our PAH-rich and PAH-poor sources. The PG QSO sample have the 
same 7.7 /xm feature strength relative to the 5.8 /xm as our PAH-poor sources, but show an 
excess of 1.6 /xm emission. 



Upcoming NICMOS observations of the sample would be able to address this directly 
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Fig. 12. — The IRAC color-color plot of our sample, excluding sources with more than one 
band at < 3 a. The colors as usual mean PAH- rich (red crosses) and PAH-poor (blue circles). 
The dashed lines indicate the wedge used in selecting AGN-dominated sources at lower- 2; 
jlacv et all 12004( 1. 
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Lastly, these results have imphcations for AGN-selection methods such as the IRAC 
color-color diagram (see Lacy et al. 2004, Stern et al. 2005, Sajina, Lacy, Scott 2005, 
Lacy et al. 2006). These rely on the power-law continuum in the IRAC band of AGN, and 
the PAH, and stellar bump features in the spectra of starbursts. In Sajina, Lacy, & Scott 
(2005), we stressed that this might fail at ^ ~ 2 due to the unknown fraction of AGN whose 
host galaxies dominate at A = 1.6/xm. Figure [12] shows the IRAC color-color plot for our 
sample. About 20% of the continuum-dominated sources fall outside of the AGN regime 
due to host galaxy emission as discussed above. This assumes that the sources undetected 
in more than 1 IRAC band and therefore not included in Figure [12] are undetected due to 
steeply falling red continua and would have fallen inside the AGN-wedge. Since our sample 
is biased toward less prominent hot dust, this is an upper limit for the AGN population as 
a whole. It does however, imply that the use of this technique to search for obscured AGN, 



although remarkably effective at lower-z (e.g. iLacy et al.ll2007l ) becomes less reliable at z>2. 



5.4. The 9.7 /xm silicate feature 



The depth of the 9 . 7 //m silicate fe ature has recently been highlighted (ILevenson et al. 



20071 : ISpoon et al.ll2007l : llmanishil 120071 ) as a tool to distinguish between sourc es dominated 



by bu ried nuclei from sources with more clumpy, extended dust distributions. ISpoon et al. 



( 120071 ) presented an intriguing diagnostic diagram of the PAH6.2 equivalent width vs. the 
9.7/im silicate feature depth. Sources with weak silicate features, range from pure starbursts 
to pure AGN in the manner of traditional PAH-over-continuum diagnostics. A distinct 
secondary sequence of increasing silicate feature depth with decreasing PAH feature strength 
is revealed. This is interpreted as an increasingly dominant buried nucleus. 

For low signal-to-noise, high-z samples such as ours, the equivalent diagram is based 
on the 7.7/im PAH feature strength, and the optical depth estimate (Figure [T3|) . The 
overlaid grids are generated via Equation 1, with an a = 2 continuum, and starburst fractions 
spanning between 3c273 (0%) and NGC7714 (100%). As we increase the optical depth the 
observed fractions remain constant, while the intrinsic fractions curve since the un-extincted 
L5.8 is higher than the observed one. We mark rg.y = 1 as separating optically thin from 
optically thick sources, which is also where extinction correction effects become strong. This 



l argely corresponds to the separation between class 1 and class 2 sources in ISpoon et al. 



(I2OO7I ) , or in other words between sources with clumpy dust (weak silicate feature) or buried 
nucleus sources (significant silicate feature). 
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Fig. 13. — The low-r, x-axis starburst fraction span the range from pure starburst 
(NGC7714) to pure typel AGN (3c273). Maintaining the same observed fraction with in- 
creasing r gives us the sohd hues. Maintaining the same intrinsic fraction with increasing r 
gives the dotted hues. The stars are the Brandl et al. (2006) starbursts. 
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Roughly 40% of the sample have rg.y < 1, indicating that a buried nucleus is not present 
or significant. Of these sources, three quarters are continuum-dominated (most likely AGN) 
sources. A quarter are strong PAH sources, but none of them are likely to be pure starbursts 
and most likely have significant AGN contribution as well. 

The bulk of our sample, ~60%, have rg.y > 1, suggesting that the buried nucleus dom- 
inates over more clumpy dust geometries. Although about 25% of these sources have sig- 
nificant PAH emission, after extinction correction ~ 90% of these highly obscured sources 
are continuum-dominated. It is important to realize that since our sensitivity to weak PAH 
emission (even using the 7.7/im feature) is worse than for the local sample used by Spoon et 
al., our diagram is missing (or are lost in the upper limits) the weakest PAH sources (class 
A in Spoon et al.). The 50% line roughly corresponds to the separation between class B and 
class C sources in Spoon et al. 

Lastly, we compare d the relatiye distr ibution of our sources along the r-direction with 



that of local ULIRGs (jSpoon et al.l 120071 ). We find that, the two are roughly consistent 
with both populations having the bulk of their members residing in the moderate optical 
depth regime (r9.7~l-3). For both local and high- 2; ULIRGs (including those with only 
lower limits on rg.y), ~25% of the pop ulation resid e s in t he most obscured sources regime 



Tgj>3). This corresponds to class 3 in lSpoon et al.l (120071 ). which contains the sources with 



evidence of extremely cold and dense media (see Appendix A). We estimate that for the local 
sample, the space density of these extremely obscured sources to be roughly ~ 6 x 10~^ Mpc~'^ 
(assuming that the ULIRG sample in Spoon et al. (2006) is not biased in its r-distribution 
compared with the IRAS IJy ULIRG sample). Given that our sample of z~2 extremely 
obscured sources is incomplete, we derive a lower limit density of ~7 x lO^^Mpc^^. These 
estimates are uncertain to within a factor of 2, but clearly suggest that the density of such 
objects was about an order of magnitude larger at 2; ~ 2 than today. 



6. Summary &; Conclusions 

We presented spectral decomposition of a sample of 48 2; ~ 2, low resolution mid-IR 
IRS spectra. The detailed analysis, separating the contribution of PAH emission, continuum 
emission and continuum extinction, has revealed several significant characteristics of z ~ 2 
ULIRGs. 

In PAPERI, we discussed that our sources have quasar-like luminosities, and frequently 
continuum-dominanted spectra suggesting an AGN is the primary power source. We find 
that 75% of our sample have continuum-dominated spectra. However, four sources of the 
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original sample of 52 were not included here because of their featureless spectra and hence 
no redshifts. If we include these 4 sources, this fraction rises to nearly 80%. About 60% 
of these show detectable PAH emission, suggestir iK these are composite sources rather than 



pure AGN. At much higher signal-to-noise ratio, ISchweitzer et al 



((20061) find -40% of PG 



QSOs have detectable PAH. We have shown that our sources have much redder mid-IR 
colors than the PG quasar sample. Lastly, in optically-selected quasars, the AGN dominates 
in the near-IR, whereas for us the host galaxy, or a significantly obscured AGN, appear to 
dominate at ~ 1 - 2 fim . The greater obscuration is also seen in the fact that ~ 60% of these 
continuum-dominated sources have rg.r > 1 (see discussion below) , whereas optically-selected 
quasars have rg 7 ~0 (and even show silicate emission). All of these suggest that the bulk of 
our sources are potentially related to the still elusive type-2 quasars. Follow-up observations 
are needed to ascertain their nature. 

About 25% of our sample have strong-PAH. In particular, for 38 spectra covering the 
rest-frame 6 /im six sources have detectable 6.2 fim PAH emission with EWq,2 ranging ~ 0.3- 
0.5 fim. Their 14 /im monochromatic luminosities are on the order of lO^^L©. ULIRGs with 
such a high infrared luminosity and such a large 6.2 /xm PAH equivalent width do not exist i n 
the local Universe. As shown by the IRS observations of local ULIRGs (lArmus et al.ll2007al ). 



a galaxy with 6.2 /im PAH equivalent width of 0.3 /im has 14 /im luminosity ~ IO^^Lq, an 
order of magnitude less than that of our z ~ 2 PAH emitting sources. Regardless of the 
exact degree of AGN contribution, it is obvious that the local PAH6.2 fim vs. Ljr relation 
(Armus et al. 2007 in prep.) does not hold at 2~2 (or experiences strong evolution). 
These 2 ~ 2 strong-PAH sources likely all have significant AGN contribution to their mid- 
IR emission, as evidenced by excess hot dust emission. However, their bolometric emission 
is likely still dominated by the starburst (by analogy with sources such as NGC6240, and 
the observed hints of colder dust at longer wavelengths). Far-IR observations are needed 
to disentangle the star-formation and AGN contribution to their infrared luminosity (see 
Sajina et al. 2007 in prep.). If these sources are indeed starburst-dominated, translating 
their continuum luminosities into bolometric luminosities (see PAPERI) implies SFR on the 
order of 1000M (7)/yr (this could be lowered significantly by including an obscured AGN, or 
top-heavy IMF (IBaugh et al.1120051)'). Such extr e me values have be en claimed before for sub- 
mm galaxies (SMGs; see e.g. iBlain et al.l 120021 : iPope et al.ll2006al ). which have comparable 
number density at 2~2 as our sources (PAPERI). 

The bulk of these luminous sources show significant obscuration with optical depths 
Tg 7 > 1. We propose a new diagnostic diagram, apphcable to low signal-to- noise ratio high-z 
data, where we plot the optical depth against the traditional PAH s trength indicator o f 
relative starburst/AGN activity. Following Levenson et al. 2006 and ISpoon et al.l (120071 ). 
about 60% of our sample are consistent with being dominated by buried nuclei rather than the 



-47- 



clumpy dust distribution indicating (but not exclusively) extended star-formation activity. 
About 50% of our strong- PAH sources fall in this category. However, we find some 2;~2, 
LiR ~ 10^^ L0 strong-PAH sources that are potentially consistent with extended, clumpy 
star- formation (e.g. MIPS2 89). Recent observations suggest this may be the case with some 



SMGs as well (jivisoni 120071 ). The definitive answer to this question will have to await the 
greater s patial resolutio n in the sub-mm regime of the Atacama Large Millimeter Array 
(ALMA; IWootten|[2006l l 



Among our sources, the distribution of the absorption depth at 9.7 fim is similar to 
what observed among local ULIRGs, with (T9.7) ~1.5; however, our sources on average are 
a factor of 5-10 more luminous (I/g-iooo/im) than local ULIRGs, and close to the bolometric 
luminosity of the brightest, optically selected quasars. The relative starburst-AGN contri- 
bution to the mid-IR emission of these sources is difficult to assess, especially for the most 
extremely obscured sources (r9.7~3), which represent ~25% of our sample . The strong 



similarity between the mid-IR spectrum of the nascent starburst NGC1377 (IRoussel et al. 



20061 ) and the deep silicate feature ULIRG IRAS08572+3915 suggest that even in sources 
with weak or absent PAH, a buried nucleus can be powered by either AGN or a young 
starburst (or more likely both). High resolution mid-IR spectra with Neon line ratios, deep 
X-ray spectra, L-band (3-4/im) spectroscopic diagnostics could provide useful additional 
information , but uncertainties due to high obscuration would remain (see e.g. discussion 



on Arp220, ISpoon et al.ll2004l ). If this deeply obscured phase is one of the ULIRG/quasar 
evolutionary stages, this phase is unlikely to last longer than ~ 10^ - 10''yrs, or until the 
supernovae- or AGN-driven winds punch holes through the cacoon. Establishing the fre- 
quency of such sources, compared with the clearly starburst-like and AGN-like ULIRGs, any 
possible evolution with redshift and luminosity would be the first steps to test whether this 
is a necessary evolutionary stage that all ULIRGs pass, or an alternate pathway depending 
on the properties of the merging components. From our sample, we estimate that 2; ~ 2 ex- 
tremely obscured sources, have a space density of at least ~ 7 x 10~^ Mpc~'^ ster^^ or about 
an order of magnitude grater than locally. 

Our study hints at some of the specific physical proc esses behind the stroiig evolution of 



the in frared luminosity function between z ~ 2 and today flLe Floc'h et al.ll2005l : ICaputi et al. 



20071 ). Although our biased selection precludes a robust comparison, we find that z^2 
ULIRGs have roughly consistent spectral properties (PAH strength, mid-IR slope, and SiO 
feature depth) as local ULIRGs, but are an order of magnitude more luminous. Thus the 
implied luminosity evolution does not follow local trends with luminosity, which would for 
example preclude the existence of strong-PAH sources with Lir~1O^^L0. The fraction of 
extremely obscured sources remains constant, but their space density increases by about an 
order of magnitude suggesting that number evolution plays a role as well as luminosity evo- 
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lution. Lastly, whether our sources are scaled-up local ULIRGs, or driven by fundamentally 
different modes of star-formation and/or black hole accretion remains an open question. 
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Fig. 14. — The absorption profile of MIPS464 after removing the best-fit continuum. We 
overlay the Galactic Center (GC) extinction curve. There is noticeable excess at ~6-7/im. 
This is most likely due to water ice and HAG absorption features. For comparison we overlay 
the 5 - 8 /im optical depth profile of NGG4418 scaled to the 9.7 /im depth of MIPS464 (which 
is ~85% that of NGG4418). Despite the low signal-to- noise ratio of our data, the agreement 
between the two suggests significant ice absorption in MIPS464. The dip at A ~ 8 is due to 
potential weak PAH component in MIPS464. The dip at A ~ 10 fim is likely due to H2 S(3) 
emission. 
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A. Water ice and crystalline silicate features in z ~ 2 ULIRGs 



About 25% of our sou rces show ex t reme silicate absorption ( t>3) consist e nt wi th the 
deeply obscured sources in ISpoon et al.l (120021 ) . or class 3A, 3B in ISpoon et al.l (120071 ). The 
high-resolution data available for the local sources discussed in Spoon et al. (2006) show that 
this class of sources frequently show additional absorption features consistent with extremely 
cold, dense media. These include in particular ices and crystalline silicates. Here we discuss 
evidence that these might be present in the most deeply obscured of our 2; ~ 2 sources. 



As mentioned in § 14. a number of our sources show ~ 6 - 7 /im absorption in addition 
to the extinction curve. Most of these are within the 1 a uncertainty of the data. The 
source in which this is most prominent is the highly-obscured MIPS464 (rg.y ~4.9, z = 1.85). 
Figure [H] shows the optical depth profile of M1PS464 compared with the Galactic Center 
extinction curve. In addition, we overlay the 5-8yum optical depth profile of NGC4418 
showing the broad 6.2 /im water ice absorption feature as well as the narrower HAC absorp- 
tion features at 6.9 /xm and 7.3 /xm. The NGC4418 profile has been scaled down by 15% 
to match the MIPS464 optical depth at 9.7/im. Although the signal-to-noise of our data 
does not allow us the detailed feature fitting that can be done for NGC4418, the good agree- 
ment between the two optical depth profiles suggests that M1PS464 is consistent with having 
water ice absorption as well. Recent studies show that water ice abs orption is obser ved in 
most local ULIRGs with high silicate feature optical depths (t9 7>3) (llmanishil 120071 ). The 
infe rred 6.0 /im abs orption depth of M1PS464 (see Figure fT^ is about twice the average of 
the llmanishil (120071 ) sample meaning that this feature is likely to be at < 1 cr for most of our 
highly obscured sources. 

We also examined our highest signal-to-noise spectra for the presence of crystalline 
substructure on the otherwise smooth 10 /im amorphous silicate feature. In the general ISM 
of our galaxy the timescale for transformation of crystalline silicates into amorphous silicates 
is estimated to be short compared w ith the timescale for injection. This is in agreement with 
the strict upper l imit of < 0.01 for the crystalline fraction of the ISM in the line of sight to our 
Galactic Center (IKemper et al.ll2004j ). Recently. ISpoon et al.l (120061 ) found the first evidence 
for crystalline silicates in 12 local ULIRGs. Surprisingly, they found a crystalline fraction 
of ~0.1, or about an order of magnitude higher than in the Milky Way disk. Figure [T5] 
shows the silicate feature profile of MIPS15880 (r9.7~3.9, z = 1.64), which has the highest 
SNR of our sample. The observed profile is remarkably similar to that of IRAS08572+3915, 
including excess absorption at 11 /im which for IRAS08572+3915 has been identified with 
a crystalline silicate feature. Given the errorbars, this opacity bump is only ~ 2 cr. The 
11.3 /im PAH feature has not been removed, in order to emphasize the significance of this 
potential excess opacity. 
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Figures [T^EiTS] show that the signal-to-noise ratio of our data does not allow detection 
of these weak features. However, the good agreement between their observed absorption 
profiles and those of two of the most obscured local sources: NGC4418 and IRAS08572+3915 
suggests similar dust properties and source geometry. 
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Fig. 15. — The silicate feature profile of MIPS15880 after removing the best-fit continuum. 
Extinction curves at best only approximate the observed profile. The IRAS08572 spectrum 
with its best-fit continuum removed, but scaled to the MIPS15880 optical depth. The two 
show very similar Si feature profiles. The discrepancy at A < 8 is due to weak PAH component 
in MIPS15880, whereas IRAS08572 does not require any PAH. The PAH have not been 
removed from the MIPS15880 profile in order to highlight the excess opacity at 11 /im even 
before the 11.3 /im PAH feature is removed. 
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Table 1. Broadband data and best- fit continuum parameters with la uncertainties. 



Source 


z 


c Ft! 

'53. 6 [MJyJ 


o r T 1 


c Ft! 


F T 1 

58 [/iJy] 


Ft! 


0. 


T9.7 


MIPS42 


1.95 


8± 3 


30± 1 


103±15 


680-1-15 

\J\J\J -1— At_J 


12.8±1.0 


1.5±0.1 


0.7±0.1 


MIPS78 


2.65 


6± 3 


39± 4 


57-1-24 


268±19 


1.5±1.3 


1.3±0.1 


3.4±0.5 


MIPSllO 


1.00 


44± 3 


43± 1 


72±18 


200-1-1 5 


9.1it3.0 


2.1±0.2 


0.5±0.5 




1.00 


57it 6 


54± 4 


55±21 




13.1±1.5 


2.3±0.1 


0.9±0.8 


MIPS180 


2.47 


23± 3 


16± 4 


14± 8 


13±17 


5.8±1.5 


1.8±0.4 


2.2±0.5 


MIPS227 


1.87 


72± 3 


82± 3 


lllilO 


113±19 


10.0±2.4 


2.9±0.2 


0.2±0.1 


MIPS27Q 


1.23 


51± 8 


52± 7 


84±15 


127±16 


9.8±1.4 


2.0±0.2 


O.OiO.l 


MIPS283 


0.94 


117± 2 


lllit 3 


147-J-23 


131±10 


10.2±2.1 


1.4±0.1 


1.5±0.9 


MIPS289 


1.86 


66± 3 


70± 4 


106±15 


93±20 


8.1±1.9 


1.9±0.3 


0.7±0.9 


MIPS429 


2.15 


18± 6 


6± 4 


4±24 


0±20 


3.5±1.5 


2.7±0.2 


>7.3 


MIPS464 


1.85 


15± 6 


19± 4 


38±16 


57±22 


8.9itl.5 


2.5±0.3 


4.9±1.6 


MIPS506 


2.52 


18± 3 


16± 1 


26±18 


13±19 


7.5±1.2 


3.6±0.2 


>6.7 


MIPS8034 


0.95 


56± 5 


64± 7 


130±23 


216±13 


6.5±1.7 


2.3±0.4 


0.3±0.3 


MIPS8184 


0.99 


104± 2 


82± 4 


82±21 


140±21 


19.7±1.4 


2.0i0.2 


0.6i0.4 


MIPS8196 


2.60 


81± 6 


61± 6 


74±19 


149±27 


5.7±1.2 


1.2±0.1 


1.3±0.4 


MIPS8207 


0.84 


117± 5 


81± 9 


94±29 


91±18 


17.9±1.6 


1.9±0.2 


0.9±0.5 


MIPS8242 


2.45 


50± 3 


43± 3 


48± 8 


88±17 


5. Oil. 5 


2.8±0.4 


0.9±0.5 


MIPS8245 


2.70 


0± 3 


9± 4 


5±15 


0±14 


1.3±1.6 


1.9±0.1 


2.8±0.7 


MIPS8268 


0.80 


95± 2 


87± 3 


94± 7 


218±14 


3.9±1.8 


1.6±0.1 


0.8±0.4 


MIPS8327 


2.48 


12± 3 


7± 4 


0±21 


46±25 


5.5±1.9 


1.9±0.1 


2.4±0.7 


MIPS8342 


1.57 


39± 6 


57± 4 


53±12 


119±19 


10.3±1.8 


2.1±0.2 


0.2±0.4 


MIPS8493 


1.80 


39± 6 


48± 4 


62±15 


24ibl4 


4.8±1.4 


1.4±0.7 


>4.7 


MIPS15840 


2.30 


18± 2 


25± 1 


62±18 


197±13 


2. Oil. 6 


1.6±0.2 


0.1±0.2 


MIPS 15880 


1.64 


63± 5 


64± 4 


58±16 


118±22 


19.5±2.1 


2.6±0.1 


3.9±0.4 


MIPS15q28 


1.52 


39± 3 


51± 4 


58±21 


125±26 


12.9±2.5 


2.6±0.2 


O.OdrO.l 


MIPS15949 


2.15 


27± 3 


30± 3 


33± 9 


89±11 


7.6±1.6 


2.3±0.2 


O.OiO.l 


MIPS 15958 


1.97 


21± 3 


43± 3 


81±24 


181±17 


8.0±2.8 


1.6±0.3 


0.5±0.4 


MIPS15977 


1.85 


57± 3 


61± 4 


69±16 


130±22 


15. Oil. 


2.2±0.2 


O.liO.l 


MIPS16030 


0.98 


84± 3 


72± 1 


94itl6 


132it20 


2.3il.6 


1.8i0.2 


0.4±0.4 


iVillr o -Luuoy 




1 7-1- '\ 


19i 4 






a 94-9 A 


9 1 -j-fi 1 


9 74-0 S 


MIPS16080 


2.04 


27± 3 


31± 7 


38±24 


71±24 


5.5±1.7 


1.7±0.1 


2.1±0.5 


MIPS16095 


1.81 


21± 2 


21± 4 


24±24 


84±14 


11.1±2.2 


2.2±0.2 


0.3±0.1 


MIPS16113 


1.90 


23± 3 


16± 3 


55±15 


35±12 


0.9±1.5 


1.4±0.4 


3.3±1.3 


MIPS16122 


1.97 


17± 3 


34± 9 


31±15 


101±15 


3.1±1.2 


1.9±0.4 


1.8±0.6 


MIPS16144 


2.13 


68± 3 


84± 3 


108±24 


105±19 


2.2il.2 


2.1±0.4 


2.0il.4 


MIPS22204 


2.08 


20± 9 


39± 6 


132±15 


478± 8 


13.2±1.1 


1.4±0.1 


1.6±0.2 


MIPS22277 


1.77 


69± 3 


67± 4 


142±15 


254± 9 


17.5±1.3 


1.9±0.2 


1.4±0.3 


MIPS22303 


2.34 


9± 3 


10± 4 


0±18 


103±13 


7.2±2.2 


1.7±0.1 


2.9±0.6 


MIPS22404 


0.61 


126± 5 


93± 3 


128±16 


117±22 


44.3±3.0 


2.2±0.1 


1.5±1.4 


MIPS22467 


0.80 


59± 9 


40± 4 


31±15 


142±25 


9.0±1.8 


1.8±0.1 


0.9±0.5 


MIPS22482 


1.84 


39± 2 


48± 3 


55±18 


93±26 


12.1±1.1 


1.9±0.2 


1.2±0.4 


MIPS22530 


1.96 


42± 5 


52± 4 


43±10 


81±12 


3.4±1.9 


2.7±0.8 


>5.2 


MIPS22554 


0.82 


99± 


79± 3 


67±13 


62±15 


11.3±1.5 


1.6±0.2 


>4.2 


MIPS22558 


3.20 


15± 3 


13± 7 


0±16 


60±16 


4.7±1.1 


2.2±0.2 


>5.6 


MIPS22600 


0.86 


93± 2 


58d= 7 


46±13 


88±20 


10.1±1.4 


1.7±0.2 


2.2±3.3 


MIPS22651 


1.73 


92± 3 


109± 4 


94±15 


99±10 


7.0±1.0 


2.0±0.3 


0.7±0.5 
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Table 1 — Continued 



Source 


z 


S3. 6 [/"Jy] 


S4.B [/iJy] 




[MJy] 


570 [mjy] 


a 


T9.7 


MIPS22661 


1.75 


,36± 2 


30± 1 


12± 8 


76±19 


2.3±2.4 


2.3±0.2 


O.OiO.l 


MIPS22699 


2.59 


9± 5 


18± 4 


17±19 


105±13 


2.5±2.3 


1.2±0.2 


1.2±0.7 
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Table 2. Luminosities with la uncertainties. Upper limits are 3a. 



Source log(Li.6/L0) log(LB.8/L0) log{Li4/LQ) log(L3o/L0) \og{L p ah ,7 .7 / Lq) 



MIPS42 









12, 


.51±0. 


.02 


12 


.68±0. 


.03 





10.60±0.20 


MIPS78 




_ 




12, 


.66±0, 


.02 


12, 


.76±0. 


,04 




<10.56 


MIPSllO 


10, 


.81±0, 


10 


11, 


.21±0, 


.07 


11, 


.64±0. 


,02 


11.84±0.52 


9.73±0.19 


MIPS133 


10 


.81±0, 


14 


10, 


.97±0, 


.05 


11, 


.45±0. 


,03 


11.71±0.27 


< 9.70 


MIPS180 


11, 


.lOitO. 


,70 


12. 


.32±0. 


.03 


12, 


.57±0. 


.07 





10.92±0.10 


MIPS227 


11, 


.53±0. 


,10 


11. 


.68±0. 


.06 


12 


.39±0. 


.02 





10.26±0.18 


MIPS279 


10, 


.89±0. 


.35 


10. 


.19±0. 


.08 


11, 


.57±0. 


.02 


11.92±0.34 


9.65±0.20 


MIPS283 


11, 


.12±0. 


,03 


10, 


.85±0, 


.08 


11, 


.01±0. 


.04 


11.48±0.72 


10.38±0.03 


MIPS289 


11, 


.47±0. 


.15 


11, 


.47±0, 


.23 


11, 


.80±0, 


.22 





11.14±0.05 


MIPS429 









11, 


.66±0, 


.09 


12 


.22±0, 


.03 





10.63±0.09 


MIPS464 








11, 


.53±0, 


.06 


12 


.16±0, 


.03 





10.09±0.19 


MIPS506 


11, 


.17±0, 


.35 


11. 


.82±0, 


.20 


12, 


.74±0. 


.01 





11.11±0.13 


MIPS8034 


10 


.75±0, 


.20 


11, 


.21±0, 


.08 


11 


.66±0. 


.02 


11.58±0.48 


< 9.59 


MIPS8184 


10, 


.88±0, 


.08 


10, 


.94±0, 


.06 


11, 


.32±0, 


.02 


11.87±0.23 


10.44±0.03 


MIPS8196 


11 


.74±0, 


.24 


12. 


.46±0, 


.02 


12, 


.53±0, 


.02 





11.09±0.08 


MIPS8207 


10, 


.89±0, 


.12 


10. 


.71±0, 


.07 


11, 


.06±0, 


.04 


11.59±0.32 


10.20±0.03 


MIPS8242 


11, 


.51±0, 


.21 


12. 


.07±0, 


.04 


12, 


.74±0. 


.04 





10.70±0.19 


MIPS8245 




— 




12. 


.23±0. 


.08 


12, 


.55±0. 


.03 


— 


10.41±0.18 


MIPS8268 


10 


.79±0. 


.05 


10. 


.82±0, 


.05 


11, 


.04±0. 


.02 


11.19±0.80 


< 9.05 


MIPS8327 









12, 


.07±0, 


.07 


12 


.40±0. 


.04 





10.62±0.21 


MIPS8342 


11, 


■ OOitO, 


.35 


11, 


.50±0, 


.08 


11, 


.93±0, 


.02 


12.28±0.36 


10.52±0.10 


MIPS8493 


11, 


.27±0, 


.27 


11. 


.38±0, 


.23 


11, 


.49±0, 


.22 





10.93±0.06 


MIPS15840 


11, 


.22±0, 


.20 


12. 


.28±0, 


.03 


12, 


.50±0, 


.04 





<10.29 


MIPS15880 


11, 


.28±0, 


.18 


11. 


.94±0, 


.03 


12, 


.51±0. 


.01 





10.60±0.09 


MIPS15928 


10, 


.95±0. 


.23 


11. 


.44±0. 


.05 


12, 


.04±0. 


.02 


12.36±0.37 


10.69±0.05 


MIPS15949 


11, 


.16±0. 


.35 


11. 


.94±0. 


.04 


12, 


.45±0. 


.03 


— 


10.71±0.08 


MIPS15958 









11, 


.95±0, 


.05 


12 


.19±0. 


.04 





10.29±0.20 


MIPS15977 


11, 


.36±0, 


.15 


11. 


.82±0, 


.05 


12, 


.28±0, 


.02 




<10.24 


MIPS16030 


10 


.98±0, 


.07 


10, 


.93±0, 


.07 


11, 


.22±0, 


.06 


11.39±0.75 


9.91±0.08 


MIPS16059 


11 


.17±0, 


57 


12, 


.05±0, 


.04 


12, 


.46±0, 


.02 




10.76±0.09 


MIPS16080 


11, 


.19±0, 


.42 


11. 


.93±0, 


.04 


12, 


.18±0. 


.02 




10.20±0.17 


MIPS16095 


10 


.85±0. 


.32 


11. 


.69±0, 


.05 


12 


.12±0. 


.02 




<10.16 


MIPS16113 








11, 


.77±0, 


.08 


11 


.88±0. 


.16 




10.80±0.08 


MIPS16122 








11, 


.81±0, 


.06 


12 


.13±0. 


.09 




10.41±0.21 


MIPS16144 


11 


.67±0. 


.10 


11, 


.59±0, 


.28 


12 


.01±0, 


.14 




11.20±0.06 


MIPS22204 








12, 


.50±0, 


.02 


12, 


.63±0. 


.01 




ll.OlitO.lO 


MIPS22277 


11 


.36±0. 


.13 


12. 


.03±0. 


.08 


12 


.36±0. 


.05 




10.44±0.19 


MIPS22303 








12. 


.32±0, 


.06 


12 


.57±0. 


.03 




<10.50 


MIPS22404 


10 


.84±0, 


.04 


10, 


.32±0, 


.07 


10 


.77±0. 


.04 


11.52±0.33 




MIPS22467 


10 


.52±0, 


.26 


10, 


.57±0, 


.05 


10 


.86±0, 


.02 


11.31±0.62 




MIPS22482 


11 


.19±0, 


.15 


11, 


.71±0, 


.08 


12 


.04±0, 


.14 




10.86±0.13 


MIPS22530 


11 


.33±0, 


.26 


11, 


.48±0, 


.25 


12 


.10±0, 


.15 




11.02±0.10 


MIPS22554 


10 


.91±0, 


.03 


10, 


.50±0, 


.09 


10, 


.71±0. 


.08 


11.29±0.57 


9.96±0.04 


MIPS22558 








12. 


.49±0, 


.02 


12 


.88±0. 


.01 




10.96±0.12 


MIPS22600 


10 


.83±0. 


.10 


10. 


.59±0. 


.10 


10 


.83±0. 


.08 


11.34±0.53 


10.04±0.05 


MIPS22651 


11 


.50±0. 


.10 


11. 


.50±0, 


.15 


11 


.86±0. 


.05 




10.77±0.07 
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Table 2 — Continued 



Source 


log(Li.6/L0) 


log(LB.8/L0) 


log(Li4/L0) 


Iog(L3o/L0) \og{LpAH,7.7/LQ) 


MIPS22661 


10.94±0.21 


11.44±0.07 


11.93±0.04 


— 10.21±0.21 


MIPS22699 




12.12±0.04 


12.17±0.08 


— <10.34 
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Table 3. PAH strength with 1 a uncertainties. Upper hmits are 2 a. 



Source 


EW6.2 
//m 


EW7.7 

/LOTl 


EW11.3 


F6.2 

10~^^ergs cm~^ 


Fr.r 

10~^^ergs cm~2 


i^ll.3 

10~^^ergs s~lcm~2 


MIPS42 


<0.08 


0.08± 0.04 


<0.04 


<6.8 


6.2±2.9 


<2.2 


MIPS78 


<0.04 


<0.06 


— 


<2.6 


<2.8 


— 


MIPSllO 


<0.26 


0.18± 0.08 


<0.08 


<6.0 


4.4±1.9 


<2.0 


MIPS133 


<0.30 


<0.30 


<0.26 


<5.4 


<5.2 


<3.4 


MIPS180 


0.08±0.04 


0.26±0.07 


— 


<2.6 


7.3±1.6 


— 


MIPS227 


<0.22 


0.17±0.08 


0.08±0.04 


<3.4 


3.2±1.3 


<2.0 


MIPS279 


— 


0.17±0.08 


0.25±0.08 


— 


2.2±1.0 


3.3±1.1 


MIPS283 


— 


2.66±0.33 


1.85±0.37 


— 


22.8±1.6 


9.9±1.6 


MIPS289 


0.48±0.14 


2.81±0.23 


0.82±0.13 


4.4±1.4 


24.3±2.7 


5.4±1.7 


MIPS429 


<0.30 


0.41± 0.09 


<2.4 


<3.2 


5.2±1.0 


<2.6 


MIPS464 


<0.12 


0.51± 0.10 


<0.32 


<1.4 


6.9±1.2 


<1.2 


MIPS506 


0.30±0.14 


0.87±0.36 


— 


3.4±1.7 


10.8±3.3 


— 


MIPS8034 


<0.26 


<0.22 


0.23±0.10 


<5.8 


<5.4 


5.8±2.5 


MIPS8184 


— 


1.83±0.11 


1.08±0.13 


— 


22.8±1.3 


11.7±1.3 


MIPS8196 


<0.06 


0.31±0.08 


— 


<2.0 


9.5±1.7 


— 


MIPS8207 


— 


2.13±0.21 


0.77±0.13 


— 


20.0±1.5 


6.4±1.1 


MIPS8242 


<0.12 


0.21±0.10 


— 


<2.4 


4.5±2.0 


— 


MIPS8245 


<0.16 


0.09±0.04 


— 


<3.4 


1.8±0.8 


— 


MIPS8268 


— 


<0.12 


<0.36 


— 


<1.6 


<3.6 


MIPS8327 


<0.14 


<0.30 


— 


<2.4 


<3.6 


— 


MIPS8342 


<0.16 


0.58± 0.14 


0.11± 0.04 


<2.2 


8.4±2.0 


1.6± 0.5 


MIPS8493 


0.32±0.13 


2.72±0.43 


2.72±0.52 


2.4±1.1 


16.4±2.2 


4.2±1.0 


MIPS15840 


<0.10 


<0.06 


<0.22 


<3.2 


<2.0 


<5.0 


MIPS15880 


0.06±0.03 


0.23±0.05 


<0.06 


<2.2 


8.9±1.9 


<1.0 


MIPS15928 


0.33±0.07 


0.89±0.11 


0.20±0.03 


4.7±1.0 


14.2±1.7 


4.0±0.6 


MIPS15949 


<0.14 


0.31± 0.06 


<0.12 


<2.8 


6.4±1.2 


<2.6 


MIPS15958 


<0.14 


0.14± 0.06 


<0.12 


<3.0 


3.0±1.4 


<1.8 


MIPS15977 


0.18±0.07 


<0.14 


0.06±0.03 


3.7±1.4 


<3.0 


1.3±0.7 


MIPS16030 


— 


0.60±0.12 


0.58±0.11 


— 


7.0±1.3 


5.6±1.1 


MIPS16059 


<0.14 


0.29±0.07 




<2.6 


5.3±1.2 




MIPS16080 


<0.14 


O.llit 0.05 


<0.20 


<2.8 


2.2± 0.9 


<1.8 


MIPS16095 


<0.12 


<0.16 


<0.06 


<2.4 


<3.2 


<1.0 


MIPS16113 


<0.16 


0.79± 0.16 


<0.38 


<2.6 


10.6± 2.1 


<1.6 


MIPS16122 


<0.16 


0.25± 0.11 


0.19± 0.09 


<2.8 


3.9± 1.9 


1.7± 0.8 


MIPS16144 


<0.32 


2.50± 0.26 


0.50± 0.24 


<2.8 


20.2± 2.8 


<2.6 


MIPS22204 


<0.08 


0.22± 0.05 


<0.10 


<7.8 


13.6± 3.1 


<2.8 


MIPS22277 


<0.12 


0.15± 0.07 


<0.08 


<4.2 


5.5± 2.4 


<2.8 


MIPS22303 


<0.12 


<0.10 




<4.6 


<3.2 




MIPS22404 






1.68±0.16 






12.7±1.2 


MIPS22467 






<0.26 






<2.2 


MIPS22482 


<0.16 


0.90± 0.15 


<0.12 


<2.6 


13.0±3.8 


<3.2 


MIPS22530 


0.39±0.15 


1.85±0.32 


0.91±0.29 


3.4±1.3 


16.2±3.8 


3.4±1.5 


MIPS22554 




2.09±0.31 


4.59±1.13 




12.0±1.1 


9.1±1.6 


MIPS22558 


<0.04 


0.16±0.08 




<1.2 


4.3±1.2 




MIPS22600 




1.94±0.32 


2.38±0.63 




13.0±1.4 


8.6±1.7 
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Table 3 — Continued 



Source 


EW6.2 


EW7.7 


EW11.3 


F6.2 


Ft. 7 


Fll.3 




//m 


//m 


IJ,m 


10~^^ergs cm~2 


10~^^ergs cm~2 


10~^^ergs s~lcm~2 


MIPS22651 


0.28±0.10 


1.12±0.18 


0.21±0.06 


3.2±1.2 


12.4±2.0 


2.0±0.5 


MIPS22661 


<0.28 


0.32± 0.16 


0.17± 0.06 


<2.6 


3.1±1.5 


1.8±0.6 


MIPS22699 


<0.12 


<0.12 




<2.0 


<1.8 
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